Evaluation of oyster exposure experiments

knitr::opts_chunk$set(echo = TRUE, warning = F, message = F, fig.path = 'figs/', dev.args = list(family = 'serif'))

library(tidyverse)
library(patchwork)
library(ggbeeswarm)
library(psycho)
library(lmerTest)
library(hrbrthemes)
library(kableExtra)
library(shiny)

source('R/funcs.R')

data(allexp)
data(bymods)
data(wtmods)

Number of folds

Raw data

yvar <- 'Number of folds'

toplo <- allexp %>% 
  filter(var %in% yvar)

p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) + 
  geom_boxplot(outlier.shape = NA, alpha = 0.1) + 
  geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) + 
  theme_ipsum_rc() + 
  facet_grid(trt~ .) + 
  scale_colour_ipsum() + 
  scale_fill_ipsum() + 
  theme(legend.title = element_blank()) + 
  labs(x = 'Exposure week', y = yvar)

p

Models by week

Olympia oyster

bymods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Pacific oyster

bymods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Models with week

Olympia oyster

dat <- wtmods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar)

dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Number of folds vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.74 0.18 4 32.45 11.48 0.00
week 0.01 0.01 2 30.48 0.38 0.69
dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Number of folds vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.63 0.16 4 25.69 9.92 0.00
week 0.01 0.01 2 24.16 0.42 0.66
trt:week 0.07 0.01 8 25.46 0.58 0.79
dat %>% 
  pull(plomod)

Pacific oyster

dat <- wtmods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar)
         
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Number of folds vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.66 0.17 4 38.14 8.31 0.00
week 0.01 0.00 2 36.80 0.21 0.81
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Number of folds vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.73 0.18 4 31.25 9.19 0.00
week 0.01 0.01 2 30.19 0.36 0.70
trt:week 0.23 0.03 8 31.05 1.42 0.23
dat %>% 
  pull(plomod)

Number of irregularites

Raw data

yvar <- 'Number of irregularities'

toplo <- allexp %>% 
  filter(var %in% yvar)

p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) + 
  geom_boxplot(outlier.shape = NA, alpha = 0.1) + 
  geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) + 
  theme_ipsum_rc() + 
  facet_grid(trt~ .) + 
  scale_colour_ipsum() + 
  scale_fill_ipsum() + 
  theme(legend.title = element_blank()) + 
  labs(x = 'Exposure week', y = yvar)

p

Models by week

Olympia oyster

bymods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Pacific oyster

bymods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Models with week

Olympia oyster

dat <- wtmods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar)

dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Number of irregularities vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.31 0.08 4 37.07 1.68 0.18
week 0.06 0.03 2 31.63 0.64 0.54
dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Number of irregularities vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.38 0.10 4 26.01 2.24 0.09
week 0.07 0.03 2 23.62 0.80 0.46
trt:week 0.53 0.07 8 26.22 1.54 0.19
dat %>% 
  pull(plomod)

Pacific oyster

dat <- wtmods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar)
         
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Number of irregularities vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.45 0.11 4 38.34 4.00 0.01
week 0.02 0.01 2 36.72 0.32 0.73
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Number of irregularities vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.70 0.18 4 27.51 6.08 0.00
week 0.03 0.02 2 26.70 0.53 0.60
trt:week 0.78 0.10 8 27.19 3.40 0.01
dat %>% 
  pull(plomod)

Total observations (dissolution)

Raw data

yvar <- 'Total observations'

toplo <- allexp %>% 
  filter(var %in% yvar)

p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) + 
  geom_boxplot(outlier.shape = NA, alpha = 0.1) + 
  geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) + 
  theme_ipsum_rc() + 
  facet_grid(trt~ .) + 
  scale_colour_ipsum() + 
  scale_fill_ipsum() + 
  theme(legend.title = element_blank()) + 
  labs(x = 'Exposure week', y = yvar)

p

Models by week

Olympia oyster

bymods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Pacific oyster

bymods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Models with week

Olympia oyster

dat <- wtmods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar)

dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Total observations vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 1.16 0.29 4 39.12 1.19 0.33
week 0.14 0.07 2 35.48 0.30 0.75
dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Total observations vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 1.42 0.36 4 31.93 1.48 0.23
week 0.16 0.08 2 29.77 0.33 0.72
trt:week 2.37 0.30 8 31.83 1.23 0.31
dat %>% 
  pull(plomod)

Pacific oyster

dat <- wtmods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar)
         
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Total observations vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 1.55 0.39 4 36.54 1.16 0.34
week 2.40 1.20 2 34.90 3.61 0.04
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Total observations vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 1.87 0.47 4 28.70 1.40 0.26
week 2.71 1.36 2 27.63 4.07 0.03
trt:week 4.51 0.56 8 28.56 1.69 0.14
dat %>% 
  pull(plomod)

Dissolution score

Raw data

yvar <- 'Dissolution score (0-3)'

toplo <- allexp %>% 
  filter(var %in% yvar)

p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) + 
  geom_boxplot(outlier.shape = NA, alpha = 0.1) + 
  geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) + 
  theme_ipsum_rc() + 
  facet_grid(trt~ .) + 
  scale_colour_ipsum() + 
  scale_fill_ipsum() + 
  theme(legend.title = element_blank()) + 
  labs(x = 'Exposure week', y = yvar)

p

Models by week

Olympia oyster

bymods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Pacific oyster

bymods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Models with week

Olympia oyster

dat <- wtmods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar)

dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Dissolution score (0-3) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 3.90 0.97 4 39.04 9.20 0.00
week 0.07 0.03 2 36.14 0.32 0.73
dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Dissolution score (0-3) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 4.78 1.19 4 28.14 10.96 0.00
week 0.05 0.02 2 26.69 0.22 0.81
trt:week 1.37 0.17 8 27.81 1.58 0.18
dat %>% 
  pull(plomod)

Pacific oyster

dat <- wtmods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar)
         
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Dissolution score (0-3) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 3.24 0.81 4 37.38 10.71 0.0
week 0.14 0.07 2 36.11 0.95 0.4
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Dissolution score (0-3) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 3.30 0.82 4 27.73 10.76 0.00
week 0.15 0.07 2 26.72 0.97 0.39
trt:week 0.49 0.06 8 27.56 0.80 0.61
dat %>% 
  pull(plomod)

Shell weight

Raw data

yvar <- 'Shell weight (g)'

toplo <- allexp %>% 
  filter(var %in% yvar)

p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) + 
  geom_boxplot(outlier.shape = NA, alpha = 0.1) + 
  geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) + 
  theme_ipsum_rc() + 
  facet_grid(trt~ .) + 
  scale_colour_ipsum() + 
  scale_fill_ipsum() + 
  theme(legend.title = element_blank()) + 
  labs(x = 'Exposure week', y = yvar)

p

Models by week

Olympia oyster

bymods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Pacific oyster

bymods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Models with week

Olympia oyster

dat <- wtmods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar)

dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Shell weight (g) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.04 0.01 4 61.89 0.30 0.88
week 0.10 0.05 2 60.58 1.36 0.26
dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Shell weight (g) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.05 0.01 4 51.80 0.32 0.86
week 0.10 0.05 2 52.45 1.29 0.28
trt:week 0.26 0.03 8 53.41 0.86 0.56
dat %>% 
  pull(plomod)

Pacific oyster

dat <- wtmods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar)
         
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Shell weight (g) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.05 0.01 4 55.71 1.29 0.29
week 0.03 0.01 2 53.03 1.29 0.28
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Shell weight (g) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.06 0.01 4 44.24 1.45 0.23
week 0.03 0.01 2 43.73 1.35 0.27
trt:week 0.11 0.01 8 44.50 1.41 0.22
dat %>% 
  pull(plomod)

Tissue weight

Raw data

yvar <- 'Tissue weight (g)'

toplo <- allexp %>% 
  filter(var %in% yvar)

p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) + 
  geom_boxplot(outlier.shape = NA, alpha = 0.1) + 
  geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) + 
  theme_ipsum_rc() + 
  facet_grid(trt~ .) + 
  scale_colour_ipsum() + 
  scale_fill_ipsum() + 
  theme(legend.title = element_blank()) + 
  labs(x = 'Exposure week', y = yvar)

p

Models by week

Olympia oyster

bymods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Pacific oyster

bymods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Models with week

Olympia oyster

dat <- wtmods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar)

dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Tissue weight (g) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0 0 4 62.14 0.56 0.69
week 0 0 2 60.84 1.40 0.26
dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Tissue weight (g) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.00 0 4 52.47 0.60 0.66
week 0.00 0 2 53.14 1.36 0.27
trt:week 0.01 0 8 54.11 0.62 0.76
dat %>% 
  pull(plomod)

Pacific oyster

dat <- wtmods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar)
         
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Tissue weight (g) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0 0 4 53.56 0.74 0.57
week 0 0 2 50.79 0.88 0.42
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Tissue weight (g) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.00 0 4 42.28 0.92 0.46
week 0.00 0 2 41.75 0.96 0.39
trt:week 0.01 0 8 42.55 1.28 0.28
dat %>% 
  pull(plomod)

Whole weight

Raw data

yvar <- 'Whole weight (g)'

toplo <- allexp %>% 
  filter(var %in% yvar)

p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) + 
  geom_boxplot(outlier.shape = NA, alpha = 0.1) + 
  geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) + 
  theme_ipsum_rc() + 
  facet_grid(trt~ .) + 
  scale_colour_ipsum() + 
  scale_fill_ipsum() + 
  theme(legend.title = element_blank()) + 
  labs(x = 'Exposure week', y = yvar)

p

Models by week

Olympia oyster

bymods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Pacific oyster

bymods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Models with week

Olympia oyster

dat <- wtmods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar)

dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Whole weight (g) vs week + treatment
Df Deviance Resid. Df Resid. Dev
NULL NA NA 305 23.08
trt 4 0.61 301 22.47
week 2 1.05 299 21.42
dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Whole weight (g) vs week x treatment
Df Deviance Resid. Df Resid. Dev
NULL NA NA 305 23.08
trt 4 0.61 301 22.47
week 2 1.05 299 21.42
trt:week 8 1.03 291 20.39
dat %>% 
  pull(plomod)

Pacific oyster

dat <- wtmods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar)
         
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Whole weight (g) vs week + treatment
Df Deviance Resid. Df Resid. Dev
NULL NA NA 290 7.68
trt 4 0.24 286 7.44
week 2 0.15 284 7.28
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Whole weight (g) vs week x treatment
Df Deviance Resid. Df Resid. Dev
NULL NA NA 290 7.68
trt 4 0.24 286 7.44
week 2 0.15 284 7.28
trt:week 8 1.07 276 6.22
dat %>% 
  pull(plomod)

Delta length

Raw data

yvar <- 'Delta length (cm)'

toplo <- allexp %>% 
  filter(var %in% yvar)

p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) + 
  geom_boxplot(outlier.shape = NA, alpha = 0.1) + 
  geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) + 
  theme_ipsum_rc() + 
  facet_grid(trt~ .) + 
  scale_colour_ipsum() + 
  scale_fill_ipsum() + 
  theme(legend.title = element_blank()) + 
  labs(x = 'Exposure week', y = yvar)

p

Models by week

Olympia oyster

bymods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Pacific oyster

bymods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Models with week

Olympia oyster

dat <- wtmods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar)

dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Delta length (cm) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.07 0.02 4 45.09 2.28 0.08
week 0.05 0.05 1 45.08 5.73 0.02
dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Delta length (cm) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.07 0.02 4 40.94 2.15 0.09
week 0.04 0.04 1 40.95 5.41 0.03
trt:week 0.01 0.00 4 40.94 0.34 0.85
dat %>% 
  pull(plomod)

Pacific oyster

dat <- wtmods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar)
         
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Delta length (cm) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.06 0.01 4 44.26 1.67 0.17
week 0.00 0.00 1 44.26 0.22 0.64
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Delta length (cm) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.06 0.01 4 40.24 1.66 0.18
week 0.00 0.00 1 40.24 0.21 0.65
trt:week 0.03 0.01 4 40.24 0.90 0.47
dat %>% 
  pull(plomod)

Delta width

Raw data

yvar <- 'Delta width (cm)'

toplo <- allexp %>% 
  filter(var %in% yvar)

p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) + 
  geom_boxplot(outlier.shape = NA, alpha = 0.1) + 
  geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) + 
  theme_ipsum_rc() + 
  facet_grid(trt~ .) + 
  scale_colour_ipsum() + 
  scale_fill_ipsum() + 
  theme(legend.title = element_blank()) + 
  labs(x = 'Exposure week', y = yvar)

p

Models by week

Olympia oyster

bymods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Pacific oyster

bymods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Models with week

Olympia oyster

dat <- wtmods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar)

dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Delta width (cm) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.06 0.02 4 44.07 1.36 0.26
week 0.04 0.04 1 44.07 3.36 0.07
dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Delta width (cm) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.06 0.01 4 40.07 1.25 0.30
week 0.04 0.04 1 40.07 3.08 0.09
trt:week 0.00 0.00 4 40.07 0.09 0.99
dat %>% 
  pull(plomod)

Pacific oyster

dat <- wtmods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar)
         
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Delta width (cm) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.03 0.01 4 44.19 0.50 0.74
week 0.00 0.00 1 44.19 0.06 0.81
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Delta width (cm) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.03 0.01 4 40.2 0.49 0.74
week 0.00 0.00 1 40.2 0.06 0.81
trt:week 0.05 0.01 4 40.2 0.89 0.48
dat %>% 
  pull(plomod)

Length rate

Raw data

yvar <- 'Length rate (cm/d)'

toplo <- allexp %>% 
  filter(var %in% yvar)

p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) + 
  geom_boxplot(outlier.shape = NA, alpha = 0.1) + 
  geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) + 
  theme_ipsum_rc() + 
  facet_grid(trt~ .) + 
  scale_colour_ipsum() + 
  scale_fill_ipsum() + 
  theme(legend.title = element_blank()) + 
  labs(x = 'Exposure week', y = yvar)

p

Models by week

Olympia oyster

bymods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Pacific oyster

bymods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Models with week

Olympia oyster

dat <- wtmods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar)

dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Length rate (cm/d) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0 0 4 45.41 2.19 0.08
week 0 0 1 45.41 1.22 0.27
dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Length rate (cm/d) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0 0 4 41.53 2.16 0.09
week 0 0 1 41.54 1.20 0.28
trt:week 0 0 4 41.53 0.83 0.51
dat %>% 
  pull(plomod)

Pacific oyster

dat <- wtmods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar)
         
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Length rate (cm/d) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0 0 4 44.22 1.85 0.14
week 0 0 1 44.22 4.38 0.04
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Length rate (cm/d) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0 0 4 40.22 1.91 0.13
week 0 0 1 40.22 4.50 0.04
trt:week 0 0 4 40.22 1.30 0.29
dat %>% 
  pull(plomod)

Width rate

Raw data

yvar <- 'Width rate (cm/d)'

toplo <- allexp %>% 
  filter(var %in% yvar)

p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) + 
  geom_boxplot(outlier.shape = NA, alpha = 0.1) + 
  geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) + 
  theme_ipsum_rc() + 
  facet_grid(trt~ .) + 
  scale_colour_ipsum() + 
  scale_fill_ipsum() + 
  theme(legend.title = element_blank()) + 
  labs(x = 'Exposure week', y = yvar)

p

Models by week

Olympia oyster

bymods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Pacific oyster

bymods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Models with week

Olympia oyster

dat <- wtmods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar)

dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Width rate (cm/d) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0 0 4 44.22 1.18 0.33
week 0 0 1 44.22 0.96 0.33
dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Width rate (cm/d) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0 0 4 40.21 1.10 0.37
week 0 0 1 40.21 0.90 0.35
trt:week 0 0 4 40.21 0.29 0.88
dat %>% 
  pull(plomod)

Pacific oyster

dat <- wtmods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar)
         
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Width rate (cm/d) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0 0 4 44.21 0.30 0.88
week 0 0 1 44.22 10.89 0.00
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Width rate (cm/d) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0 0 4 40.2 0.29 0.88
week 0 0 1 40.2 10.46 0.00
trt:week 0 0 4 40.2 0.56 0.69
dat %>% 
  pull(plomod)

Final length

Raw data

yvar <- 'Final length (cm)'

toplo <- allexp %>% 
  filter(var %in% yvar)

p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) + 
  geom_boxplot(outlier.shape = NA, alpha = 0.1) + 
  geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) + 
  theme_ipsum_rc() + 
  facet_grid(trt~ .) + 
  scale_colour_ipsum() + 
  scale_fill_ipsum() + 
  theme(legend.title = element_blank()) + 
  labs(x = 'Exposure week', y = yvar)

p

Models by week

Olympia oyster

bymods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Pacific oyster

bymods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Models with week

Olympia oyster

dat <- wtmods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar)

dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Final length (cm) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.16 0.04 4 52.98 0.45 0.77
week 0.24 0.12 2 52.95 1.35 0.27
dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Final length (cm) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.20 0.05 4 44.90 0.54 0.71
week 0.26 0.13 2 44.91 1.44 0.25
trt:week 1.05 0.13 8 44.91 1.44 0.21
dat %>% 
  pull(plomod)

Pacific oyster

dat <- wtmods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar)
         
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Final length (cm) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.20 0.05 4 52.96 0.56 0.70
week 0.06 0.03 2 52.91 0.33 0.72
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Final length (cm) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.14 0.03 4 44.86 0.38 0.82
week 0.06 0.03 2 44.89 0.34 0.72
trt:week 0.73 0.09 8 44.89 1.02 0.43
dat %>% 
  pull(plomod)

Final width

Raw data

yvar <- 'Final width (cm)'

toplo <- allexp %>% 
  filter(var %in% yvar)

p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) + 
  geom_boxplot(outlier.shape = NA, alpha = 0.1) + 
  geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) + 
  theme_ipsum_rc() + 
  facet_grid(trt~ .) + 
  scale_colour_ipsum() + 
  scale_fill_ipsum() + 
  theme(legend.title = element_blank()) + 
  labs(x = 'Exposure week', y = yvar)

p

Models by week

Olympia oyster

bymods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Pacific oyster

bymods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Models with week

Olympia oyster

dat <- wtmods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar)

dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Final width (cm) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.18 0.05 4 53.03 0.67 0.61
week 0.23 0.12 2 52.99 1.69 0.19
dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Final width (cm) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.17 0.04 4 44.97 0.61 0.66
week 0.23 0.12 2 44.99 1.72 0.19
trt:week 0.60 0.08 8 44.99 1.11 0.38
dat %>% 
  pull(plomod)

Pacific oyster

dat <- wtmods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar)
         
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Final width (cm) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.06 0.02 4 53.00 0.31 0.87
week 0.11 0.06 2 52.93 1.07 0.35
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Final width (cm) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 0.12 0.03 4 44.85 0.56 0.70
week 0.11 0.05 2 44.89 1.05 0.36
trt:week 0.38 0.05 8 44.89 0.91 0.52
dat %>% 
  pull(plomod)

Respiration

Raw data

yvar <- 'Respiration (umol/hr/g)'

toplo <- allexp %>% 
  filter(var %in% yvar)

p <- ggplot(toplo, aes(x = factor(week), y = val, colour = species, fill = species)) + 
  geom_boxplot(outlier.shape = NA, alpha = 0.1) + 
  geom_point(aes(group = species), position = position_jitterdodge(jitter.width = 0.1), size = 2, alpha = 0.6) + 
  theme_ipsum_rc() + 
  facet_grid(trt~ .) + 
  scale_colour_ipsum() + 
  scale_fill_ipsum() + 
  theme(legend.title = element_blank()) + 
  labs(x = 'Exposure week', y = yvar)

p

Models by week

Olympia oyster

bymods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Pacific oyster

bymods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar) %>% 
  pull(plomod)

Models with week

Olympia oyster

dat <- wtmods %>% 
  filter(species == 'Olympia') %>% 
  filter(var == yvar)

dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Respiration (umol/hr/g) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 16.68 4.17 4 51.04 0.34 0.85
week 9.67 4.83 2 55.99 0.40 0.67
dat %>% 
  filter(species == 'Olympia') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Respiration (umol/hr/g) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 13.93 3.48 4 44.51 0.29 0.88
week 10.23 5.12 2 46.96 0.42 0.66
trt:week 64.95 8.12 8 47.14 0.67 0.72
dat %>% 
  pull(plomod)

Pacific oyster

dat <- wtmods %>% 
  filter(species == 'Pacific') %>% 
  filter(var == yvar)
         
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomodnoint) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week + treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Respiration (umol/hr/g) vs week + treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 40.25 10.06 4 53.53 0.73 0.58
week 11.31 5.66 2 52.35 0.41 0.67
dat %>% 
  filter(species == 'Pacific') %>% 
  pull(anomod) %>% 
  .[[1]] %>% 
  knitr::kable(format = 'html', digits = 2, caption = paste0(yvar, ' vs week x treatment')) %>%
  kable_styling(full_width = T, font_size = 14) %>% 
  HTML
Respiration (umol/hr/g) vs week x treatment
Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
trt 37.54 9.39 4 45.50 0.69 0.61
week 9.84 4.92 2 44.45 0.36 0.70
trt:week 54.96 6.87 8 45.48 0.50 0.85
dat %>% 
  pull(plomod)